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1. Introduction 



The hydrogen atom in constant homogeneous electric and magnetic fields still remains to be 
an object of theoretical investigations. A good example is the recent work [|IJ, where a recurrent 
nonperturbative method is developed for building the exact wave function of hydrogen atom in 
magnetic field in the form of convergent double series. More wide discussion of the problem is 
contained in the review 0. 

The famous technical trouble, namely the inability to separate the variables, only stimulate 
application of new investigation methods, including perturbative ones. The moment method ||, 
first used for perturbative treatment of the anharmonic oscillator, is not related with variables 
separation. The recent application of this method to the Zeeman effect problem [|]] allowed to check 
the behavior of high orders asymptotic of the perturbation series. The moment method in the 
form similar to that used in ||], was developed independently of Ader's work || by Fernandez and 
Castro [H . Then it was applied to hydrogen atom placed in parallel electric and magnetic fields || 
and later the Zeeman effect problem was considered for four sequences of hydrogen atom states @. 

It seems to be even more important to apply it to hydrogen atom in crossed electric (S) and 
magnetic (7i) fields because only initial terms of expansion in powers of £ and 7i were considered 
for this case up to now ||, [], [H], 12]. As will be shown here the moment method allows to 



compute high enough orders of this expansion. 



The high orders asymptotic can be obtained with the help of the imaginary time method |L3, [14 



lq| . This asymptotic is determined by the contribution of an extreme subbarrier classical trajectory 
into the atom ionization probability |16|, |17|]. A pair of extreme paths replaces this trajectory at 
some value of the ratio of external fields 7 = 7i/£. The dependence of far perturbation series terms 
on 7 reflects this change of extreme trajectory and should be especially sharp for perpendicular 
external fields. We study here just this case. 

2. Recurrent evaluation of perturbation series 

Consider the ground state of hydrogen atom, placed in perpendicular electric and magnetic 
fields. These fields are supposed to be constant and homogeneous. We restrict ourselves with 
nonrelativistic approximation and neglect the spin of electron. From the very beginning we take 
measures to simplify the numerical computations and to achieve high enough order of perturbation 
theory. For this aim we consider 7 as a fixed parameter, replacing the double expansion in external 
fields by the single-variable series 

00 00 

V = e = y,e^, (i) 

where the wave function corrections \k) and hyper-polarizabilities depend on 7. We introduce 
also circular coordinates 

x± — x ± iy , 
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then all further relations will have real coefficients. In these coordinates the hamiltonian of our 
problem is: 

H = H + £Hi + £ 2 H 2 ; = ~A - - ; 

2 r 

^=x + Il z = -{x + + x„) + 1[x + —-x„—} ; (2) 

2 2 

H 2 = — (x 2 + y 2 ) = —x + x_ . 
The wave function correction of the order k satisfies the differential equation 

[k/2] 

(H - E )\k) = -H^k - 1) - H 2 \k - 2) + E 2j\k ~ 2j) . (3) 

3=1 

Just as in other problems where the moment method was used |3|, M it is not difficult to transform 
equation (3) into algebraic relation between moments of the order k 

PL/3 = {Q\r a - a -Px%xt\k) . (4) 

A recurrence relation results 
I 
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o-a-P) (a + « + + l)P^ + 2apP*_ 2>a _ 1>f3 _ 1 - = R k ^ , (5) 



where 



1 'Y 2 

2 1 

[k/2] 



Paaj3 — [P(r+l,a+l,P + P<r+l,a,P+l + l{ a P)Paap] + o ^a+2, 



- V E L P k ~ 2j 

i=i 

The right-hand side of eq. (5) and hyper-polarizability depend only on the moments of preceding 
orders. As usual in the moment method ||, the orthogonality condition is accepted 

(0|fc) = <f Ol * — ^ ^o,o = V- (6) 
An expression for hyper-polarizability arises from eq. (5) at o = a = (3 = and even k 

Ei = i(P&3 + Pto'x) + jP2 k xi ■ (7) 

The closed system of recurrence relations (5) - (7) allows to achieve, at least in principle, an arbitrary 
high order of perturbation theory. The sequence of operations is similar (also somewhat simpler) 
to that, used in the work Q to compute Zeeman's shift of a non-degenerate state. At every order k 
only moments P k a p from the sector a > a-\- (3 — 1, a>0,/3>0 are necessary. They are evaluated 
by successively increasing of a, a and (3 values with the help of eq. (5). 

We have obtained hyper-polarizabilities in perpendicular fields up to 80th order, see table 1. 
This order is large enough to compare the dependence of these coefficients on 7, see fig. 1, with the 
predictions, following from quasiclassical considerations. One can see from fig. 1, that the function 
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fki'j) = ha.(\E^-\/k\) has two features. It has a minimum at 7 ~ 3.4 and a sequence of singular 
points to the right of this value. Besides, the function £^(7) changes it's sign at every singular 
point of f k (7). 

As follows from table 1, at not very large 7 values all coefficients have negative sign, as in 
the case of Stark effect. In intermediate region of 7 values the sequence of E^ signs is irregular and 
for sufficiently large 7's the series has normal Zeeman's sequence of signs (— 1) J+1 . 

3. High orders asymptotic 

As is well known ||16|| , a dispersion relation connects asymptotic of high orders coefficients Efc 
with the ionization probability of the atom i.e. with the penetrability of the potential barrier. This 
relation arises as a consequence of the fact, that the energy eigenvalue E = E (£ 2 ) — fT(£ 2 ) has 
essential singularity at £ 2 = and a cut along £ 2 > semiaxis. (And similarly E{7i 2 ) has essential 
singularity at Ti = and a cut TL 2 < 0.) 

To evaluate the ionization probability V the imaginary time method was previously developed |13 



14], The leading term of the asymptotic E^ of E^ coefficients at k — > 00 is determined by the 



classical subbarrier path with extremal value of the abbreviated action. Time takes complex values 
during this subbarier motion. There are two kinds of complex classical trajectories. Like in the Stark 
effect case, the ionization may be caused by electric field, at stabilizing influence of the magnetic 
field. The path of this kind creates the asymptotic 

£^(7) ~ k\ a fc ( 7 ) , k is even , (8) 



at not very large magnetic field, for 7 below some critical value j c . According to [Fj| 7 C = 3.54 
for perpendicular external fields. And it is possible to cross the barrier also at Ti 2 < 0, like in the 
Zeeman effect problem. Subbarrier trajectories of this kind are responsible for the form of £^"(7) 
in the opposite case 7 > j c . This change of asymptotic explains the origin of the left minimum in 
fig. 1. 



Having in mind to get estimate for the function 0(7), entering E£, we apply the results of [18|, [19 1 
and write here some necessary expressions for the special case of perpendicular external fields. More 
general considerations related to arbitrary £ and 7i mutual orientation are contained in the work |T8|. 



The time of subbarrier motion satisfies the equation [19 



r 2 - (rcthr - l) 2 = 7 2 (9) 

which has a set of solutions r n = irnr + r^. The minimal value of the imaginary part of the subbarrier 
action is provided by tq for 7 < 7 C and by a pair of solutions t±i for 7 > j c . In the region 7 < 7 C 
the energy half-width is 



T(£ ) = -^exp 



2s(t) 



3£ 



3r 
2^ 



9(1) = 7^ h 2 - • (10) 
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The dispersion relation in £ 2 then leads to 



E. 



2.7 



1 °° 

= ~2^I 



T(z)dz 



2.7 



where 



a(i) 



2<7(7) 



(11) 



(12) 



The last equality is valid also in the region 7 > j c , where #(7) and 0(7) are complex functions. At 
7 < 7 C the resulting approximate expressions for 0(7) are 



0(7) 



30 



7 



0(7) 



71 
" 2100 

4 7 3 



7* + 



( 7 2 _ i)2 (1 _ 2e -7 2 -i) 
And in the region 7 > 7 C another representation works 

'8tt 2 



7< 1 ; 
7 > 1 • 



K7)l = - 

7T 



1 "-2 + 

7 2 



+ 3 



1 



(13) 
(14) 

(15) 



On the other hand in the limit of large k a simple relation appropriate for numerical evaluation 
holds: 

\Ek\ 



In (0(7) I = — In- 
dk 



k\ 



(16) 



Evaluating 0(7) above j c , we used smoothed function £^-(7), with the nodes vicinities excluded. A 
comparison of this way numerically obtained function 0(7) with expressions (13) - (15) is presented 
in fig. 2. 

Now we turn our attention to the region 7 > 7 C . Two solutions of eq. (9) T\ and r_i lead to 
complex conjugate values of (7(7). Substituting approximate T\ value into second expression (10), 
it is possible to get the phase of the function 0(7): 

arg (a) = -arg (g) = -- + 0(7), 0(7) = - - 

2 7 

Finally the sign-alternating asymptotic arises: 



3 7 3 



+ 0(l/ 7 5 ). 



E. 



2.7 



2|5( 7 )|(2j)!|a| 2 ^ 1 



cos 



7T 



(2j + l)(-- + «( 7 ))+/5(7) 



(17) 



~ (2j)\ \a\ 2j+1 sin [(2j + l)a( 7 ) + /3( 7 )] , j » 1 . 

Here ^(7) = ar(? (B) is the phase of the preexponential factor. Its relative contribution to the total 
phase falls like 1/j. 

When the order of perturbation 2j is fixed and 7 increases, expression (17) changes its sign at 
every point where the argument of the sinus turns to zero. This could explain the singular points 
in fig. 1 in the language of asymptotic. But rather lengthy calculations are required to establish 
detailed quantitative correspondence between asymptotic (17) and exact E^ coefficients, including 
nodes vicinities. Simple approximate expression for a (7) is not enough for this aim. 
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4. Discussion 



For the general case of the ground state energy expansion in powers of crossed external fields, 
the term of the forth power was known long enough II . 

E=- 1 - + ±E^; = + ; (18) 

E{A) = _3|5 - 4 + 1|^ 2 + 10 m , _ 53^ (ig) 

The value of is confirmed for perpendicular fields by the work [12] and for parallel fields - 
by m^EH. The E£ coefficient, computed by means of recurrence relations (5) - (7) exactly agree 
with (19). But we have noticed numerical difference between our coefficient E§ and corresponding 



quantity from the work [ 12]] . Therefore the sixth order of perturbation theory was analyzed in 
details. 

The magneto-electric susceptibilities, i.e. coefficients of the double series in powers of external 
fields, can be easily obtained from hyper-polarizabilities £^(7). Thus, in the sixth order, taking 
into account that Stark's and Zeeman's coefficients are fixed, it is enough to choose four different 7 
values and to solve the system of four linear equations. The following representation results 

= J2l6-2j,2j(K/£) 2j ■ (20) 
(The last identity introduces notation of fl2|| .) Using linear relation between expansions (1) and (18) 



and the known magneto-electric susceptibilities in parallel fields p0 |, it is easy to obtain another 
term of series (18): 

,~ 2512779 - 254955 - 4 -o 133199 ^ r -- l2 
512 512 256 1 J 

_49195 g, - 4 _ 255557 H - 2|7?£? + 5581^ (21) 
1536 6912 L J 4608 v ; 

Some next terms of series (18) can be obtained in the same way. Expressions (20) and (21) are 

convenient to check term by term the sixth order correction. As follows from |12| 



1610197 ± pi 2417015 



27648 !J 1536 

while the results of our computation are 



953869 , , 521353 

794. = and 7 49 = . (23) 

12A 13824 142 512 v ; 



All other corresponding coefficients of [12| and of present work coincide. We carried out additional 



independent calculation by means of the method from the work and get 



, fl = _953869 

m 13824 ' V ; 
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see Appendix. Note, that || contains complete correction of the sixth power in external fields 
for the case of parallel fields and only a part of it for the case of perpendicular fields. These 
"celebrated" sixth order terms result as a by-product of forth-order calculations in the work ||. 
The agreement between high-order hyper-polarizabilities Efc and their asymptotic E^ presents 
additional confirmation of correctness of recurrence relations (5) - (7). 

5. Concluding remarks 

The considered above problem demonstrates once more the high efficiency and convenience of 
the moment method. The obtained recurrence relations have allowed to advance up to 80th order of 
perturbation theory. Besides the unusual "oscillations" of hyper-polarizabilities as a function of the 
ratio of external fields were noticed. The high orders asymptotical behavior was analyzed as well. 
Basic parameters of this asymptotic exactly agree with those, previously obtained on the ground of 
quasiclassical approximation with the help of imaginary time method. 
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Table 1. Hyper-polarizabilities Et of the hydrogen ground state. 
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Appendix 

Extending the described in the work calculations we obtained, by the same method, the 
ground state energy correction, which is proportional to Q = 7i 2 [HS] 2 . The perturbation in || 
includes Stark's term V s = £r, paramagnetic V p = \'HL and diamagnetic Vp = ^[Hr\ 2 terms. 
The entire perturbation is inhomogeneous, therefore terms of the sixth power in external fields are 
presented in corrections of the fourth, fifth and sixth orders in V. 

151347 



r(4) 



-Q + ... 



(Al) 



r (5) 



r (6) 



2047 

(2\(V - e (1) )|2) - 2e (2) (2|l) - e (3) (l|l) , (A2) 

r (2) ((2|2) + (1|3» - 2e (3) <2|l> - e (4) (l|l) . (A3) 

In the following an abbreviated mnemonic notation will be used, reflecting the origin of each term 
and the powers of entering this term external fields. This notation helps to omit all not essential 
terms. Operation signs are encircled in the abbreviated notation. In the first order in V 



(3\(V-eW)\2) 



|1> ^ { ai (r)(Sr) + -(a 2 (r)H 2 + a 3 (r)[HS] 2 )}\0) B{V £ ®Vd}\0) 



(A4) 
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The next correction |2) contains 



Vp^|O)0^ 1 (r)([^>)|O) ) 



VeV D \V>) e±{(Sr)(b 4 (r)[Hr\ 2 + b 5 (r)H 2 ) + b e (r)H 2 ([HS}r)}\0) . (A5) 

The polynomials a«(r) and bi(r) are given in the article 0. Abbreviated notation allows to verify 
that in each of the right-hand sides of Eqs. (A2) and (A3) only the first matrix element yields 
contribution ~ Q. 

e& = A + B 1 + B 2 + C+ ... , 

dots stand for all omitted terms. 



AS(0\{VpVe)V D {V P V e )\0) , 
Bie{0\(V e V D )V P (VpVe)\0) , 
B 2 e{0\(VpV e )Vp(V e V D )\0) , B 2 = B l 
Ce-eW(0\(V P V £ )(VpV £ )\0) . 



(A6) 



(A7) 



Corresponding explicit expressions are 

A = ^m{r){[n^ 2 [Hff |0> = , 
B i = ^\bi(r){h(r)[Hf>\ 2 + (b 5 (r) + & 6 (r))ft 2 }([7?£]f) 2 |0> 

1 18432 9673 

Only one term of the third correction to the wave function is essential - that of the lowest power 
in external fields: 

|3> = ipi(r)([W[7^]TO|0> + ...e(^%e..O|0> • (AS) 
The differential equation for ci(r) || is satisfied by the polynomial: 

ci(r) = -^( 450 + 225r + 62r2 + 6r3 ) • (^ 9 ) 

As a consequence we get 

e^B(0\(V 2 V s )Vp(VpV £ )\0} + ... (A10) 

and 

£ (6) = _l(o| &1 ( r ) Cl ( r )([7Y[7i£]]r) 2 |0) + ... . (Alt)') 
16 

The total energy correction of the desired form is 

AE = -^^ln 2 [HS] 2 + ... . (All) 
6912 1 J K J 



One should not forget also the "isotropic" contribution to the energy correction, originating from 

49195 



1536 



-H 4 £ 2 . 
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Figure captions 

\E ± \ 

Fig. 1. Functions fk{l) = In k resulting from the recurrently computed hyper- 

k\ 

polarizabilities. 

Fig. 2. Parameter 0(7) of the perturbation series asymptotic. 
The solid line follows from the quasiclassical estimate at 7 <C 1, see eq. (13); the same estimate for 
7 ^> 1 is presented by dashed lines, see eqs. (14) and (15). Numerically obtained values are denoted 
by stars. 
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